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Abstract. Recent work in percolation has led to exact solutions for the site and 
bond critical thresholds of many new lattices. Here we show how these results can be 
extended to other classes of graphs, significantly increasing the number and variety of 
solved problems. Any graph that can be decomposed into a certain arrangement of 
triangles, which we call self-dual, gives a class of lattices whose percolation thresholds 
can be found exactly by a recently introduced triangle-triangle transformation. We 
use this method to generalize Wierman's solution of the bow-tie lattice to yield several 
new solutions. We also give another example of a self-dual arrangement of triangles 
that leads to a further class of solvable problems. There are certainly many more such 
classes. 
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1. Introduction 

The simple geometrical model of percolation on regular lattices has been studied 
extensively since it was first proposed nearly fifty years ago p. It continues to find broad 
applications to such diverse problems as understanding conductive materials jSHHj, the 
fractality of coastlines [I], networks El Ej, turbulence (HI Ej, colloids [TQJ , and the 
spin quantum Hall transition ^T] . Along with the Ising model of a ferromagnet and the 
related Potts model, percolation serves as a paradigm of statistical mechanics. 

In the usual bond percolation model each edge (bond) of an infinite graph or lattice 
is defined to be open with probability p (and closed with probability 1 —p) independently 
of all other bonds. If p is near most bonds are closed, but there are small and 
sparsely distributed open clusters. With p close to 1, the landscape is dominated by 
an infinite open cluster. Between these two regimes, there lies a sharply defined critical 
threshold, p c , which is a characteristic of each lattice, at which the infinite cluster 
makes its first appearance. Finding this critical point is central to many theoretical 
and computational studies of percolation. Yet, for many years, exact thresholds were 
known for only a relatively small number of two-dimensional lattices (square, triangular, 
honeycomb, bow-tie...) found through duality and a technique called the star-triangle 
transformation [T2*] I13j. 

Recently, it has been recognized that the star-triangle transformation can be 
generalized for systems with correlated bonds jTU H3j or for systems composed of 
triangular cells more complicated than a simple star or triangle [10]. Using this 
generalization, site (in the site problem we consider vertices occupied with probability 
p) and bond thresholds have been found for many new lattices, including the "martini" 
(Figure E]), "martini- A" and "martini-B" lattices [Tl] llfij. Recently, Wu has shown 
that the transition point of the Potts model (which includes both the Ising model and 
percolation as special cases) can also be found explicitly for these three lattices - the 
first new exact transition points for the Potts model found in many years. 

In this report we briefly discuss these new results, and show how they can be 
extended further to systems where the triangular cells do not necessarily lie in a 
triangular lattice, but in any self-dual arrangement (in a sense defined below). In this 
way we provide a derivation of the threshold for Wierman's bow-tie lattice ^3] and 
generalizations of it, as well as other new lattices. 

2. Triangle-triangle transformation 

We consider any system that can be decomposed into distinct triangles, like those shown 
in Fig. ^ For lattices of this type that are also invariant under a certain transformation, 
which we call the triangle-triangle transformation, it is possible to determine exactly the 
bond threshold. The triangle-triangle transformation is illustrated in Figure ^a) where 
we superpose a reversed (hatched) triangle over each (shaded) triangle of the original 
lattice (Fig. E^c)). If the lattice formed by the hatched triangles is the same as the 
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Figure 1. The triangle-triangle transformation is shown in (a). The lattices in (b) 
and (d) are invariant under this transformation, as shown in (c) and (e). 




a) b) 

Figure 2. a) The kagome lattice. As clear from b) this system is not self-dual under 
the triangle-triangle transformation. 

original, we say that the lattice is self-dual under the triangle-triangle transformation. 
The shaded triangle (the generator) can represent any network of bonds, correlated 
bonds, and sites, as long as they are contained between three vertices. This means 
that every self-dual arrangement of triangles gives an infinite family of solvable lattices. 
Figure El shows an example, the kagome arrangement, that can be decomposed into 
triangles but is not self-dual under the triangle-triangle transformation. Although some 
special lattices fall into this arrangement, such as the double bond honeycomb lattice, 
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Figure 3. Some examples of lattices created by using different generators in the 
arrangement of Fig. E^b). a ) The martini lattice, b) The martini- A lattice, c) the 
martini-B lattice. 



and can be solved exactly, this is not true for every choice of generator. For example, 
choosing a simple triangle with uncorrelated bonds in Figure |2f a), gives the kagome 
lattice. This lattice's bond threshold cannot be found by the following method and, in 
fact, the problem remains unsolved [TH] . 

To find the general condition for criticality, consider any self-dual arrangement with 
vertices A, B, C on the original triangles, and A', B', C on the dual triangles as in Fig. 
^ To be at the critical threshold, we require two things (in exact analogy with the usual 
star-triangle transformation): (1) that the connection probabilities between the vertices 
of A, B, C and those of the vertices of A', B', C are the same, and (2) that A, B, C 
and A', B', C are dual to each other. The former implies that the overall connection 
probabilities between points on the lattice and corresponding ones on the (identical) 
dual lattice are the same, and the latter implies that this equi-probability point must 
be the critical point (because only at the critical point can a lattice and its dual have 
identical correlations). 

To state these conditions mathematically, we let P(A, B, C) equal the probability 
that A, B, and C on a triangle are connected, P(A, B, C) equal the probability that A 
and B are connected but not C, etc., and P(A, B,C) equal the probability that none 
are connected. Likewise define P'(A', B', C) etc. for the connection probabilities on the 
dual triangles. Then the condition of equality of connection probabilities is expressed 
by 



P(A,B,C) = P'(A',B',C') 



(1) 
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P{A,B,C) = P'{A',B',C') 
P(A,B,C) = P'(A',B',C') 



(2) 
(3) 



Secondly, we require that the two triangles A, B, C and A', B', C are dual to each other. 
This requires 



For example, the first says simply that if A, B, and C are connected to each other, then 
A', B', and C must all be isolated. This is because the dual triangular unit A', B', C is 
created by drawing bonds intersecting all the bonds of A, B, C, and placed down with 
complementary probilities Pi = 1 — pi. Thus, whenever there is an occupied bond on 
the original lattice, the dual-lattice bond is not occupied, and vice-versa. Therefore, if 
A, B, and C are all connected, the bonds between A', B' and C must necessarily be all 
unconnected. The same kind of argument holds for the other conditions. 

Comparing the two sets of equations (HJ)-© an d Q-®, we see that the middle 
(two-point) equations are identical. The first and third equations yield the requirement 
that the probability that all the vertices connect is the same as the probability that 
none of them connect: 



This simple general relation determines the critical point, and was given independently 
in Refs. [15] and [16]. It is only a function of the arrangement of bonds in the individual 
triangles and not dependent upon how the triangles are arranged to form a lattice, as 
long as they are arranged in an overall self-dual way as those in Fig. 1. In our previous 
work [TB| , was applied only to systems of triangular cells that were on a triangular 
array (Fig. db)). The lattices we studied, along with their generators, are shown in 
Figure El As already mentioned, the site thresholds for these lattices can also be found 
exactly. In fact, they were initially derived by a different method ^1] but were then 
shown to fall into the current framework by using generators that involve sites as well as 
bonds [in] • These bond and site thresholds are given in Tables [T] and El Here we apply 
(J7j) to any self-dual set of triangles such as those in Figs. [T](d) and|Hl The generalization 
of (JIJ) to all self-dual arrangements is immediately obvious from the arguments of Refs. 
[inillS], and was suggested by the earlier work of Wierman on the bow-tie lattice [T3*] . 

3. Martini and related lattices 

The result from which all our exact solutions will follow is the critical threshold of the 
martini lattice in which each bond of the generator is assigned a different probability, 
as shown in Figure IH^a). This results in a critical surface rather than a critical point 
and application of (JZj) implies that it is given by [El HHj 



P(A,B,C) 
P(A,B,C) 
P(A,B,C) 



p'(a',b',c') 

P'(A',B',C) 
P r (A',B r ,C) 



(4) 
(5) 
(6) 



P(A,B,C)=P(A,B,C). 



(7) 



1 - Pij9 2 r 3 - p 2 P3ri - pmr2 ~ V\V2T\r% 
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Figure 4. Some examples of generalized bow-tie lattices, 
lattice p^ ond equation 

martini 1/^^ 0.707107... (2p 2 - l)(p 4 - 3p 3 + 2p 2 + 1) = 
martini-A 0.625457... p 5 - 4p 4 + 3p 3 + 2p 2 - 1 = 

martini-B 1 /2 (2p - l)(p 2 - p - 1) = 

Table 1. Bond percolation thresholds of the lattices in Fig. [21 and their associated 
polynomials derived from equation (JHJ. 

- PiP3nr 3 - p 2 p 3 r 2 r 3 + P\P2Vz r \ r i 
+ PiP2P3nr 3 + p 1 p 2 P3r 2 r 3 + pip 2 rir 2 r 3 

+ PiP3nr 2 r 3 + p 2 p 3 r x r 2 r 3 - 2p 1 p 2 p 3 r 1 r 2 r 3 = . (8) 

If we set all bonds equal to p, we get the condition (2p 2 — l)(p 4 — 3p 3 + 2p 2 + 1) = 0, 
which has solution on [0, 1] p c = l/y/2. Setting Pi — 1 and p 2 = p 3 = ri = r 2 = r 3 = p 
gives the martini-A lattice, and setting p 2 = p 3 = 1 and p\ = r\ = r 2 = r 3 = p gives the 
martini-B lattice, as listed in Table ^ 

Taking r*j = 1, we get the result for the critical surface of the general honeycomb 
lattice [I2|, 1 — Pip 2 — P1P3 — p 2 p3 + P\p 2 P3 = , and taking Pi — 1 yields the critical 
surface of the general triangular lattice IT? . 1 — r\ — r 2 — r 3 + rir 2 r 3 = . Setting 
all probabilities equal in the last two expressions gives polynomial equations with the 
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Figure 5. The duals of the lattices shown in Fig. 0] 

lattice p^ lte equation 

martini 0.764826... p 4 - 3p 3 + 1 = 

martini-A l/y/2 « 0.707107... 2p 2 - 1 = 
martini-B 1/0 « 0.618034... p 2 + p - 1 = 

Table 2. Site percolation thresholds of the lattices in Fig. and their associated 
polynomials. <f> = (1 + V / 5)/2 = 1.618034... is the golden ratio, which makes an 
interesting appearance here. 

well-known solutions p c (honeycomb) = 1 — 2sin7r/18 and p c (triangular) = 2sin7r/18. It 
is a general fact that thresholds found via (J7J) are solutions of polynomials with integer 
coefficients. 

4. Bow-tie class of lattices 

Although equation (jHJ) was originally derived for the martini lattice, it gives the critical 
surface of any lattice in which the martini generator is arranged in a self-dual way. Here 
we apply the result to the arrangement given in Fig. ^d). Taking pj = 1 in this case 
gives the lattice of Fig. ^d) where the shaded triangle is just a simple triangular cell. 
This lattice has a double bond where the bases of abutting triangles join. However, 
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Figure 6. a) Another self-dual triangular lattice under the triangle-triangle 
transformation, b) a particular example of a) using the "A" lattice generator 



lattice 


„bond 
Pc 


equation 


a) 


0.404518 


1 - p - Qp 2 + Qp 3 - p 5 = P3] 


b) 


0.533213 


1 - p - 2p 3 - Ap A - 4p 5 + 15p 6 






+ 13p 7 - 3Qp 8 + 19p 9 + p 10 + p u 


c) 


0.672929 


1 - 2p 3 - 2p 4 - 2p 5 - 7p 6 + 18p 7 






-flip 8 - 35p 9 + 21p w - Ap 11 = 


d) 


0.625457 


1 - 2p 2 - 3p 3 + Ap 4 - p 5 = 


e) 


0.595482 


= 1 -p c {a) 


f) 


0.466787 


= l-Pc(b) 


g) 


0.327071 


= 1 -Pc(c) 


h) 


0.374543 


= l-Pc{d) 



Table 3. Thresholds of the bow-tie lattices in Figs.0]and[S] 



we can use a trick, due to Wierman j!3j . to combine these two bonds by setting the 
probability of each to be r 3 = 1 — ^/l — p. The probability of traversing the double 
bonds is now just p, which effectively combines them into one. With this replacement, 
we get Wierman's bow-tie lattice (Fig. Ufa)). Setting the single bonds, r x and r 2 , to p 
gives the condition 1 — p — 6p 2 + 6p 3 — p 5 = 0, with solution in [0, 1] p c = 0.404518... 
as found by Wierman. But we can just as well make the basic triangular cell the "A" 
generator (Fig. Efb)) which produces the lattice shown in Fig. Iljd), with a transition 
point identical to that of the regular "A" lattice (given in Table HJ) since no double bond 
arises in this case. Of course, there is an infinite variety of additional lattices that can be 
solved by this procedure and it becomes a somewhat subjective matter to decide which 
ones are most interesting; some of our choices are shown in Fig. 01 with their duals in 
Fig. Eland thresholds in Table El For (c), we use a martini generator with pi = ^Jp and 
P2 = P3 = n = ^2 = r 3 = p, while for (b) we use an "A" generator with an extra bond 
on the bottom, which results in a doubled bond that must be combined. 
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In Fig. Efa) we show another arrangement of triangles that is self-dual under 
the triangle-triangle transformation, defining another family of solvable lattices. For 
example, replacing the shaded triangles by the "A" generator, one can get an interesting 
variation on the "A" lattice with the same threshold, shown in Fig. EJ^b). Many similar 
self-dual arrangements can be constructed. 

Thus, we have shown that a great variety of lattices with exactly determined 
percolation thresholds can now be constructed. We also see that there are many 
variations of a given lattice (such as the lattices with the A-lattice threshold shown 
in Figs. 3(b), 4(d), and 6(b)) that share the same threshold, even though the bonds are 
in different arrangements. Even after decades of study, the percolation model continues 
to yield new exact results. 
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